%lambda为λ，dt为Δti，M为求到的最大正整数
function fL1=findfL(lambda,dt,M)
fL1=zeros(1,M+1);%同样，多给一列表示0的情况，后面依次顺延
for m=0:M
    fL1(m+1)=exp(-lambda*dt)*(lambda*dt)^m/factorial(m);
end